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6.1 DISCRETE TIME-FREQUENCY DISTRIBUTIONS 0 

For the purposes of digital storage and processing, any real-life signal that is not 
discrete-time and time-limited must be made so by sampling and windowing. More- 
over, if we wish to evaluate a continuous time- frequency distribution (TFD) numeri- 
cally in a finite number of operations, we must be content with computing a sampled 
form of the TFD from a finite number of samples of the signal. For such reasons, 
we need to define discrete, time-limited equivalents of continuous TFDs. This ar- 
ticle derives discrete forms of the Wigner-Ville distribution (WVD), the windowed 
WVD and the general quadratic TFD, and gives thirteen examples of discrete-time 
kernels. Thus it extends the material present in Chapters 2 and 3. 

If the signal z(t) is ideally sampled at times t — n/ f s , where n is an integer and 
f a is the sampling rate, it becomes 

OO OO 

z(t) E *(*-£) = E *(£)*(*-£) ( 6 - L1 ) 

n=—o o n =— oo 

where S (. . .) denotes the unit impulse function. We shall use a wide caret Q to 
denote a TFD that has been modified by sampling. In the lag (r) domain, we shall 
consider only ideal sampling at r = m// s , where m is an integer. This in turn will 
draw attention to the discrete time values t = n/ f s , where n is an integer; but it 
will not be necessary to specify ideal sampling in the time domain. 

6.1.1 The Discrete Wigner-Ville Distribution (DWVD) 

The WVD of a continuous-time signal z(t) has the form 

/ OO 

z{t+l)z*{t-^)e-^ T dr. (6.1.2) 

-OO 

Changing the variable of integration to 6 = r/2 , we obtain the alternative definition 

/ oo 

z(t + 6)z*{t-e)e- jA7rf0 de (6.1.3) 

-oo 

which we shall find more convenient for the purpose of conversion to discrete form. 
The product z(t+9) z*(t — 9) is the instantaneous autocorrelation function (IAF) 
in terms of t and 9. Let us consider this IAF as a function of 9 and suppose that it 
is sampled. 

Theorem 6.1.1: If W z (t, f) is modified by ideally sampling the IAF in 9 at 
9 = r/2 = m/f s , where m is an integer and f s is the sampling rate, and if the 
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modified TFD is denoted by W z (t, /), then 

^(^>5#) = 2 X z ( n jT k ) e -^ km/N (6.1.4) 

\m\<N/2 



where N is the duration (in samples) of z(n/f s ); and the time support for W z is 
the same as for z(n/f s ). 

Proof/explanation: After sampling, the integrand in Eq. (6.1.3) becomes 



OO 

z(t + e) z*(t-e) e - j4nfe J2 s ( d ~%) ( 6 . 1 . 5 ) 

m= — o o 

and the WVD becomes 

OO 

W z (t,f)= 2 J2 z(t+^)z*(t-^)e~ j47rfrn/fs . (6.1.6) 

m= — o o 



Now let t be conveniently restricted to t = n/f s where n is an integer (this represents sampling 
of the signal in time, but the sampling is not ideal). Then, if z{t) has finite duration, the time 
support of z(n/f s ) can be written in the form \n — no\ < N/2 , where no is real. Hence the time 
support of W z is 

|n zb m — no | < N/2 (6.1.7) 

where both signs must be satisfied. Because the sign can be chosen so as to increase the magnitude, 
we must have \n — no \ + \m\ < N/2 , hence |n — no| < N/2 and \m\ < N/2 . So the time support 
for W z is the same as for z(n/f s ), while the summation in Eq. (6.1.6) is restricted to \m\ < N/2 , 
giving a maximum of N terms. 1 The sampling in r makes W z (t, f) periodic in / with period f s / 2 , 
while the time-limiting gives a frequency resolution of N bins per period. So it is convenient to let 



/ = 



kfs 

2N 



where k is an integer. With these restrictions, Eq. (6.1.6) reduces to Eq. (6.1.4). 



(6.1.8) 



With a change of notation, Eq. (6.1.4) becomes 

W z [n,k]= 2 ^ z[n+m\ z*[n — m\ e ~^ 2vkm / N . (6.1.9) 

\m\<N/2 



This W z [n,k] is called the discrete WVD or DWVD. If z[n] has a duration not 
exceeding N samples, the DWVD is represented as an NxN real matrix. If the 
summand is extended periodically in m with period N (i.e. extended periodically 
in r with period 2 N/f s ), we obtain 

W z [n, k\ = 2 DFT {z[n+m\ z *[n—m]} ; mG(N) (6.1.10) 

m — >k 

where ( N ) means any set of N consecutive integers. 2 



1 N terms for odd N; N— 1 terms for even N. 

2 For even N, the periodic extension is padded with a zero term. 
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It remains to find the minimum value of f s that avoids aliasing in 
the frequency and Doppler domains. As usual, we define the instanta- 
neous autocorrelation function as K z (t, r) = z(t+|) z*(t — the spectrum as 
Z(f) = T{z{t)} , the WVD as W z (t, f) = IF {K z (t, r)} , and the ambiguity func- 
tion as A z (v, t) = T {I\ z (t, r)} . Using these notations and the familiar properties 
of the FT, it can be shown that 

W z (t, f) = [2Z(2f) e^‘] * [2Z*(2f) e~W] (6.1.11) 

A z (u,t) = [Z{y) e j ™ T ] * [Z*(-v) eT jW ] . (6.1.12) 

If Z(f) is zero outside the band |/| < B/2, then the spectrum on the right of 
Eq. (6.1.11), i.e. the WVD, is zero outside the band |/| < B/2, and the spectrum 
on the right of Eq. (6.1.12) is zero outside the band \u\ < B . Aliasing will be avoided 
if the sampling rate in r (namely f s / 2, since f s is the sampling rate in 9) is at least 
B and the sampling rate in t (namely f s ) is at least 2 B\ that is, aliasing will be 
avoided if / s < 2 B . The minimum sampling rate in r makes the WVD periodic in 
/ with period B. 

The sampling rate may be reduced in the case of an analytic signal. If z(t) is 
the analytic associate of the real signal s[t), whose spectrum S(f) is zero outside 
the band |/| < B/2, then Z(f) is zero outside the band 0 < / < B/2, so that 
the spectrum on the right of Eq. (6.1.11), i.e. the WVD, is zero outside the band 
0 < f < B/2 , and the spectrum on the right of Eq. (6.1.12) is zero outside the band 
\v\ < B/2. Aliasing will be avoided if the sampling rate in r is at least B/2 and 
the sampling rate in t is at least 5; that is, aliasing will be avoided if f s <B. The 
minimum sampling rate in r makes the WVD periodic in / with period B/2. 

Now consider W z [n, k] as a matrix with k as the column index. If the sampling 
rate in t is B , there is no zero-padding in the frequency domain for an analytic 
signal, and all the columns of W z [n,k\ are needed for the positive frequencies. If 
the sampling rate in t is 2 B, only half the columns of W z [n, k] are needed to represent 
the positive frequencies. The negative-frequency elements may be assumed to be 
zero if z[n\ is analytic, provided of course that the assumption of analyticity is not 
invalidated by short-segment effects, such as windowing of the signal in the n domain 
and/or windowing of the IAF in the m domain. The latter kind of windowing is 
discussed under the next heading. 

6.1.2 The Windowed DWVD 

Nonlinear IF laws and multiple signal components give rise to artifacts (interference 
terms, cross-terms) in the WVD [see Article 4.2]. The effect of nonlinear IF laws 
can be reduced, with a concomitant loss of frequency resolution, by windowing the 
IAF in the lag direction (i.e. in r or 9) before taking the FT. For continuous time, 
if the window is 



9(t) =g ( |) =g{0) , 



(6.1.13) 
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then the resulting TFD, denoted by Wf(t, /), is 

/ OO 

9{t) z(t + f ) z*(t - f ) e-^dr (6.1.14) 

-OO 

/ OO 

g{0) z(t + 0) z*(t - 9) e- j ^ f0 d9 (6.1.15) 

-OO 

and is called the windowed WVD or pseudo- WVD. The effects of sampling and 
time-limiting on the windowed WVD are described by the following theorem. 

Theorem 6.1.2: If Wj!(t,f) is modified by ideally sampling g(r) in 9 at 
9 = t/2 = m/f s , where m is an integer and f s is the sampling rate, and if g(r) 
is time-limited so that 

g(r) = 0 for \9\ = |§| > ^ (6.1.16) 

where M is a positive integer, and if the modified TFD is denoted by Wf(t. /), then 

WHf.^) = 2 E g{f)z{^)z\^)e-^ km / M (6.1.17) 

\m\<M/2 

and the time support for Wf is the same as for z(n/f s ). 

Proof: Apart from the limits on r, 0 and m, which lead to the substitution / = f/j , the 
explanation is similar to that of Theorem 6.1.1. ■ 



With a change of notation, Eq. (6.1.17) becomes 

Wj![n,k\= 2 ^2 g[rn]z[n-\-m]z*[n — m]e~^ 2 ' Kkrn / M . (6.1.18) 

|m|<M/2 

This Wf [n, k] is called the windowed DVWD or pseudo-DWVD. If z[n] has 
a duration not exceeding N samples, the windowed DWVD is represented as an 
N x M matrix. If the summand is extended periodically in to with period M (i.e. 
extended periodically in r with period 2 M/f s ), we obtain 

WH [n, k\ = 2 DFT {g[m\ z[n+m\ z*[n— to]} ; toG(M). (6.1.19) 

m—*k 



6.1.3 The Discrete Quadratic TFD 

In terms of continuous variables, the general quadratic TFD is 



Pz(t,f) = JF f {G{t,T)*[z{t+^) z *{t-^)]} 

POO POO 

= / G(u,t) z(t—u+^) z*{t—u—^) due~^ 2 ^ T dr 



( 6 . 1 . 20 ) 



— OO 



— OO 



( 6 . 1 . 21 ) 
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where G(t, r) is the time-lag kernel. If we define 

G(t,0) = G(t,T) (6.1.22) 

where 9 ■= t / 2 , we obtain the alternative definition 

/ OO nOO 

/ G{u,9) z(t — u+9) z*(t—u—9) e-i 4n f e dudO. (6.1.23) 

- OO d — OO 



Theorem 6.1.3: If p z (t, /) is modified by ideally sampling G(u,t) in 9 at 
9 = t/ 2 = m/ f a and in u at u = p/ f s , where m and p are integers and / s is the 
sampling rate, and if G(u , r) is time-limited so that 

G(u,t) =0 for |0| = HI > 5 ^ or |u| > ^ (6.1.24) 

where M and P are positive integers, and if the modified TFD is denoted by p z {t, /), 
then 



n (H- Ml) — 2 

> 2M ) ~ Z 






bl<§ 






-p+m \ 

.fs ' 



p—m \ —j2Tvkm/M 



f n—p—m \ 

V fs ) 



(6.1.25) 



Proof/explanation: The sampled version of G(u,0) is 

OO OO 

G(u,0) = G(u,0) J2 S ( u ~jt ) J2 <5(0- ^) (6.1.26) 



p= — oo m= — oo 

OO oo 

= G(u,0) £ £ S 2 {u-f s ,0-f s ) (6.1.27) 

p= — oo m= — oo 



where 62 ( 11 , 6 ) is the two-dimensional unit impulse function. When G(u, 6 ) is replaced by G(u, 9), 
Eq. (6.1.23) becomes 



Mt,n = 2 J2 E 



G(^ 



/s’ /, 






“j47T fm/f s 



m= — 00 p= — 00 



(6.1.28) 



The time-limiting of G(u, r) restricts the summation to |m| < M/2 and |p| < P/2 . The sampling 
in r makes p z (t,f ) periodic in / with period /s/2, while the time-limiting in r (or m) gives a 
frequency resolution of M bins per period. So it is convenient to let 



/ = 



kj s 
2M 



(6.1.29) 



where /.■ is an integer. With these restrictions, and with t = n/f s , Eq. (6.1.28) reduces to 
Eq. (6.1.25). ■ 



With a change of notation, Eq. (6.1.25) becomes 

*M] = 2 £ E G[p, m] z[n— p+m] z*[n— p— m] e i 2nkm / M (6.1.30) 

"'l< a l/'l< a 

= 2 J2 G[n,m]*(z[n+m]z*[n-m])e- j2nkm / M . (6.1.31) 

\m\<f 
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This p z [n,k\ is the generalized discrete quadratic TFD. 

If the summand is extended in m with period M , Eq. (6.1.31) becomes 

/) z [n, fe] = 2 DFT {G[n,m] j (z[n+m]z*[n-m])} ; rnG(M). (6.1.32) 



Because the time (n) support for z[n+m] z*[n—m\ is the same as for z[n], the time 
support for p z [n,k\ is the same as for G[n,m] *z[n ] . If the latter has a duration 
not exceeding N samples, then the non-zero elements of p z [n, k] may represented 
as an NxM matrix and the linear convolution in Eqs. (6.1.31) and (6.1.32) may 
be interpreted as modulo- N. If G[n,m\ is real and even in m, then the argument 
of the DFT is Hermitian in m, so that p z [n, /c] is real. 

Eq. (6.1.32) shows that the implementation of the discrete quadratic TFD in- 
volves construction of the discrete IAF, followed by convolution in n with the time- 
lag kernel, followed by discrete Fourier transformation. These steps may be simpli- 
fied by taking advantage of symmetries, as explained in [1,2]. 

It remains to determine the effect of G(t , r) on the required sampling rate. If 
T {G(t, r)} = y(f, /) and T {G(t, r)} = g(v, r), we have the familiar results 

-r— ►/ t—n; 



JF f {g(t,r) *K z (t,r)} = q(t, /) ** W z (t, f) (6.1.33) 

TS 9 ^ t,T ' ) * K z (t,i-)} = g{v,r) A z {y 7 T). (6.1.34) 

Comparing the above with Eqs. (6.1.11) and (6.1.12), we see that there is an addi- 
tional spreading of the spectrum in / but not in v. If y(t, /) is zero outside the band 

|/| < Bq / 2 , then the total bandwidth of the WVD is increased by Bq, so that the 
required sampling rate is increased by Bq in r, or 2 Bq in 9 , and the sampling rate 
in t must be increased to match. 



6. 1.3.1 Special Cases 

If G[n, m ] = S[n ] g[m ] , then Eq. (6.1.32) reduces to Eq. (6.1.19), so that the discrete 
quadratic TFD reduces to the windowed DWVD. If, in addition, g [to] = 1 (that is, 
if G[n,m] =(5[?r]), then Eqs. (6.1.19) and (6.1.32) reduce to Eq. (6.1.10), so that 
the windowed DWVD and the discrete quadratic TFD reduce to the DWVD. 

Two of the three theorems above concern the sampling of a window or kernel 
function. Theorem 6.1.1 seems to be an exception in that the entire IAF is sampled 
(which is possible only in theory, as one cannot compute a continuous IAF in prac- 
tice). However, because the WVD may be considered as a windowed WVD with 
g{r) = 1 , Theorem 6.1.1 can be restated in terms of sampling the lag window, like 
Theorem 6.1.2. 
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6. 1.3. 2 Doppler-Frequency Form 



The Doppler-frequency form of the general quadratic TFD is 

Pz(t,f) = F-^{g(v,f)*[Z(f+%)Z*(f- 1)]} (6.1.35) 

/ oo /»oo 

/ g{v,v)Z{f-V+%)ZV-V-%)drie’ 2 * vt dv (6.1.36) 

-ooJ—oo 

where Giy, /) is the Doppler-frequency kernel. If we define 

e(S,f) = 0(v,f) (6.1.37) 

where £ = v/2 , we obtain the alternative definition 

/ oo poo 

/ g&tfZif-ri+OZV-ri-Oei^dvdZ. (6.1.38) 

- oo J — oo 

As the time-lag definition of the quadratic TFD leads to Theorem 6.1.3, so the 
Doppler-frequency definition leads to the following result. 



Theorem 6.1.4: If p z (t, /) is modified by ideally sampling G{ri 1 v) in £ at 
£ = v /2 = and in 77 at rj = , where l and q are integers and N is a posi- 

tive integer and f s is a positive constant, and if G(rj , v) is band-limited so that 

G{r),v) =0 for |£| = ||| > jfy- or |?y| > Qfy- (6.1.39) 

where L and Q are positive integers, and if the modified TFD is denoted by p z (t, /), 
then 



n (n_ kfs\ _ 9 
V f. ’ 2N ) ~ z 



E ^ 

*l<# 



l9l<? 



n(lfs qfa \ ry( [fc — g+Z]/ B \ ry* ( [k- q-l] f e \ „j2nln/N 
y \2N’2N/ Zj \ 2 N I V 2 N / e 



(6.1.40) 



Proof: Parallel to the proof of Theorem 6.1.3 and the subsequent discussion. ■ 

With a change of notation, Ecp (6.1.40) becomes 

p z [n,k] = 2^ G[l,q]Z[k-q+l}Z*[k-q-l]e jMn / N (6.1.41) 
ld<f kl<§ 

= 2 ^ 0[Z, jfe] * (Z[k+l] Z*[k—l}) e j2 * ln/N . (6.1.42) 

M 

If the summand is extended in l with period L , Eq. (6.1.42) becomes 

p z [n,k} = 2IDFT{G[l,k]*(Z[k+l]Z*[k-l})y, ls(L). (6.1.43) 

The frequency support for p z [n,k] is the same as for G[l,k] *Z[k \ . If the latter 
support has a width not exceeding K frequency samples, then the non-zero elements 
of p z \n,k\ may be represented as an LxK matrix and the linear convolution in 
Eqs. (6.1.42) and (6.1.43) may be interpreted as modulo- AT. While the dimensions 
of the TFD matrix seem to differ from those in Theorem 6.1.3, the dimensions are 
upper bounds and may be matched, if desired, by zero-padding. 
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Table 6.1.1: Kernel requirements for selected properties ("Prop.”) of quadratic TFDs, in time-lag 
and Doppler-frequency domains, for general continuous and discrete TFDs, and discrete TFDs with 
Doppler-independent (Dl) and lag-independent (LI) kernels. Note that 6 = t/2 and £ = v/2. 



Prop. 


KERNEL CONSTRAINTS 


Continuous 


Discrete 


Discrete DI 
G[n,m} = 5[n\g2 [m] 
G[l, k] = G 2 [k] 


Discrete LI 
G[n, m] = gi [n\ 
Q[l,k] = Gil]S[k] 


RE: 


G(t,r) 

= G*(t, -r). 


G[n, m\ 

= G*|n, — m]. 


*(N 

II 

(M 


gi [n] is real. 


TM: 


G(f, 0) = <5(t) 


II 

cT 

sf 

O 


92 [0] = 1 


WVD only 


FM: 


G(0,f) = 5(f) 


d[0,k] = <S[fc] 


WVD only 


Gi [0] = 1 


IF: 


G(t, 0) = 5(t) ; 

J fG(uJ)df = 0. 


G[n, 0] = <5[n] ; 
Z k kg\i,k]^o. 


92 [0] = 1 ; 
E fe fcG 2 [fc] =0. 


WVD only 


TD: 


G(o,f) = 5(f) ; 
f tG(t,r) dt = 0. 


G[0 ,k] = <5[fc] ; 
E„«G[n,m] = 0. 


WVD only 


Gi [0] = 1 ; 
En n SlW = °' 


TS: 


G(t, t) = 0 
if |0|<|t|. 


G[n, m\ = 0 
if \m\ < |n|. 


Always 


WVD only 


FS: 


G{y, f) = o 
if |£|<|/|. 


G[l,k] = 0 
if |i|<|fe|. 


WVD only 


Always 



6.1.4 Desirable Properties; Kernel Constraints 

The desirable properties of continuous TFDs (defined in Section 3.1.1, p. 60 ff) are 
easily redefined for discrete TFDs. Some important examples are given below. 
Realness (RE) says simply that the TFD is real. 

The marginal conditions, which may be considered optional for signal- 
processing purposes, are the time marginal (TM) 



Y, k pMM = \M\ 2 


(6.1.44) 


and the frequency marginal (FM) 




E„ P MM = \z[k]?- 


(6.1.45) 



The IF property says that the periodic first moment of the TFD w.r.t. fre- 
quency is the instantaneous frequency. Its dual, which seems to be regarded as less 
important, is the time delay property (TD), and says that the periodic first 
moment of the TFD w.r.t. time is the time delay. 

The time support (TS) property says that if z[n] = 0 everywhere except 
m < n < n 2 , then p z [n, fc] = 0 everywhere except n\ < n < ri 2 • Similarly, the 
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Table 6.1.2: Kernels of selected TFDs in time-lag and Doppler-frequency domains. For the spectro- 
gram and windowed Levin (lu-Levin) distributions, the window w[n] is assumed to be real and even. 
W[k\ denotes the sampled spectrum of the window w(r). In the property column, an exclamation 
(!) means that the property is always satisfied, while an asterisk (*) means that the property is satisfied 
subject to non-degenerate constraints on the window and/or parameter. 



Distribution 


KERNEL 


PROPERTY 

RE TM FM IF TD TS FS 


G[n , m] 


Q\l,k] 


WVD 


<S[n] 


5[k\ 


!!!!!!! 


Levin 


t<5[n+m] 

+ ^5[n — m] 


^6[k+l] 

-| ~^S[k — l] 


!!!!!!! 


Born- Jordan 


L\am\ rect (tem)] 
** [ sine n sine m] 


[ 4|aII rect ( 4^7 )] 
** [ sine k sine l] 


!!!!!** 


Modified B 


cosh - n 
Y n cosh -2 / 3 n 




! ! ! ! 


windowed WVD 


<5[n] w[m] 


W[k\ 


* * * i 


re-Levin 


7^w[m] S[n-\-m] 
J r\w[m\ S[n — m] 


\W[k+l] 

+ \W[k-l } 


i * * i 


Zhao- Atlas-Marks 


[w[m] rect(f^)] 
** [ sine n sine m] 




* * * * 


Rihaczek 


S[n — m\ 


S[k + l] 


! ! ! ! 


te-Rihaczek 


w*[—m\ S[n — m] 


W*[k + 1\ 


* ! 


Page 


S[n — \m\] 


5<5[fc + ;] + j(5[fc — Z] 

+ ^d=m* sincl 


! ! ! ! ! 


Choi-Williams 


2|m| CX P( 4m 2 ) 

** [ sine n sine m] 


v/ ^exni _,r2 ' Tfc2 'l 
2 1 Z | 4 l 2 ) 

** [sine k sine l] 


! ! ! ! ! 


B 


/ |2m| \ P 

1 . 9 ' * sine m 

\ cosh^n / 




i 


spectrogram 


w[n-\-m] w[n — m] 


W[k+l]W[k-l\ 


i 



frequency support (FS) property says that if Z[k] = 0 everywhere except 
ki < k < Jt 2 , then p z [n, k] = 0 everywhere except fci < k < &2 • 

The class of TFDs satisfying realness, the time and frequency marginals, the IF 
property and the time and frequency support properties is called class V. 

The properties listed above are selected from Table 3.3.1 on p. 75. For each 
property, Table 3.3.1 gives necessary and sufficient conditions on the kernel of the 
general continuous quadratic TFD. To obtain the corresponding conditions for dis- 
crete TFDs, we first express the conditions entirely in the time-lag and Doppler- 
frequency domains, using the scaled variables 6 and £ where convenient, obtaining 
the “Continuous” column of Table 6.1.1. Then we sample the kernels as specified 
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in Theorems 6.1.3 and 6.1.4, obtaining the “Discrete” column of Table 6.1.1. The 
remaining columns are obtained by specialization. 

The sampling of the time-lag kernel will be free of aliasing if the kernel is first 
band-limited in t and 6 to ±/ s /2. This causes S(t) to be discretized as <5[n]. Similarly, 
the sampling of the Doppler-frequency kernel will be free of aliasing if the kernel is 
first time-limited in / and £ to ±N/f s . This causes 6(f) to be discretized as S[k], 

6.1.5 Examples 

By sampling the kernels of common continuous quadratic TFDs [see Table 3.3.2 on 
p. 76], we obtain the new Table 6.1.2, which lists the kernels for the discrete versions 
of those TFDs. The convolutions with sine functions are performed before restricting 
the variable to integer values; this requires oversampling and is computationally 
inefficient. The convolutions in n and m arise from the band-limiting of the time- 
lag kernel prior to sampling. This band-limiting of the kernel does not affect the 
result of convolving the kernel with the IAF, because the IAF is assumed to be 
similarly band-limited. Similarly, the convolutions in k and l arise from the time- 
limiting of the Doppler-Frequency kernel prior to sampling. Where no convolution 
appears in the kernel, either a sine function has been converted to a discrete delta 
function by the sampling, or a window function is assumed to provide sufficient 
filtering. 

6.1.6 Summary and Conclusions 

Ideal sampling of window or kernel functions leads to straightforward definitions of 
discrete-time forms of the WVD, the windowed WVD and other quadratic TFDs. 
Use of the analytic signal minimizes the required sampling rate. 

Further theoretical details may be found in [3-5] . Some practical computational 
issues will be examined in Article 6.5. 
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